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Abstract 

Choosing a simple class of flows, with characteristics that may be present in the 
Earth's core, we study the ability to generate a magnetic field when the flow is per- 
mitted to oscillate periodically in time. The flow characteristics are parameterised 
by .D, representing a differential rotation, M, a meridional circulation, and C, a 
roll component characterising convective rolls. The dynamo action of all solutions 
with fixed parameters (steady flows) is known from previous studies. Dynamo ac- 
tion is sensitive to these flow parameters and fails spectacularly for much of the 
parameter space where magnetic flux is concentrated into small regions, leading 
to high diffusion. In addition, steady flows generate only steady or regularly re- 
versing oscillatory fields and cannot therefore reproduce irregular geomagnetic-type 
reversal behaviour. Oscillations of the flow are introduced by varying the flow pa- 
rameters in time, defining a closed orbit in the space (D, M). When the frequency 
of the oscillation is small, the net growth rate of the magnetic field over one period 
approaches the average of the growth rates for steady flows along the orbit. At 
increased frequency time-dependence appears to smooth out flux concentrations, 
often enhancing dynamo action. Dynamo action can be impaired, however, when 
flux concentrations of opposite signs occur close together as smoothing destroys the 
flux by cancellation. It is possible to produce geomagnetic-type reversals by mak- 
ing the orbit stray into a region where the steady flows generate oscillatory fields. 
In this case, however, dynamo action was not found to be enhanced by the time- 
dependence. A novel approach is taken to solving the time-dependent eigenvalue 
problem, where by combining Floquet theory with a matrix-free Krylov-subspace 
method we avoid large memory requirements for storing the matrix required by the 
standard approach. 

Keywords: Kinematic dynamos; Time dependent stability; Geomagnetism; Floquet 
theory; Eigenvalue problems. 
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1 Introduction 



Unlike the Sun or Jupiter, the Earth's dynamo runs on a tight heat bud- 
get and may therefore rely on a significant large-scale component to the 
flow. In addition to the dominant dipole component of the geomagnetic 
field, there is some observational eviden ce to suggest the field con tains a 
persistent non-axisymmetric component ( Gubbins and Kellv . 19931 ) . Both 



may be indicative of a steady component to the underlying core flow. This 
prompted a kinematic stud y of dynamo a c tion fr o m a class of large-scal e 



prompted a Kinematic stud y ot d ynamo acti on tr o m a class ot la rge-scale 
steady candidate core- fl ows llGubbins et all fcOQOal . iGubbins etafl lEOOOhh . 



Gubbins and Gibbons! ( 2002 ) hereafter referred to as papers I— III). Steady 
flows, particularly if they contain stagnation points, tend to concentrate 
magnetic flux into small regions with large energy loss due to diffusion. Of- 
ten, increasing the flow speed to overcome diffusive loss simply results in 
more concentrated flux and faster decay. On the other hand, chaotic flows 
appear to make better dynamos, perhaps because the mixing properties of 
the flow prevents permanent flux concentration and exponential separation 
of neighbouring particle s in the flow lead to stretching of the magnetic field 



f article 
1998h . The dynamo mechanisms of complicated chaotic 



flows are difficult to understand. We are therefore motivated to first study 
the effects of simple flows that fluctuate about a steady mean. Although 
the Earth's dipole has persisted for a long time, secular variation includ- 
ing excursions of the magnetic pole may indicate that fluctuations of the 
large-scale flow are present in the Earth's core. 

Kinematic theory ignores the nonlinearity of back-reaction by the mag- 
netic field on the flow, and considers only the time evolution of the magnetic 
field b as governed by the induction equation 

d t h = R m V A (u A b) + V 2 b. (1) 

The induction equation has been non-dimensionalised with the timescale for 
magnetic diffusion, the length scale d, the radius of the sphere, and in this 
work the velocity u is normalised such that the magnetic Reynolds number 
R m is unity for a flow of unit kinetic energy. For a given steady flow the 
induction equation is linear in b and has eigenfunction solutions of the form 
b(r, 0,(f>)t) = e at ~B(r,6,(j)). Dynamo action is established if 3?(cr) > 0. This 
simple test is the major advantage of the kinematic approach. The alterna- 
tive is to integrate in time the nonlinear problem for both the velocity and 
magnetic field until one is convinced the magnetic field will not ultimately 
decay; this is expensive and the results can be uncertain. The advantage re- 
mains, however, when the flow varies in time but is periodic; Floquet theory 
giv es an eigenvalu e problem for the growth rate. 



:s an eigenvalu e problem tor tne growtn rate. 

Backus! ( lgssl ) was first to show kinematic dynamo action by a time- 



dependent flow. His dynamo employed periods of stasis while high harmonics 
in the field decayed, enabling him to establish convergence of the solution. 
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Time dependence may even lead to dynamo action when no single snapshot 
of the flow can generate magnetic field on its own. Magnetic fields can 
grow during an initial transient period under the influence of a subcritical 
steady flow, the most familiar example being the production of toroidal field 
from the action of differential rotation on poloidal field. If the induction 
equation were self-adjoint and its eigenfunctions orthogonal it would be a 
simple matter to prove that all such transients decay; it is the non-normal 
property that allows transients to grow. The initial fields which optimise 
transient growth for flows in a sphere, including one of the flows here, have 
been studied by iLivermore and Jackson! (|2004f ) . Unfortunately, if the flow 



is steady the field eventually dies away and only the slowest decaying mode 
remains. If the flow is permitted to be time dependent, however, once a 
transient field associated with the initial flow has grown, a change in the 
flow can encourage further growth. In plane-layer flow it has been shown 
that by repeatedly switching t he orientat i on of the flow it is possible to take 



advantage of these transients (G og et all 119991 ). and to find dynamo action 



where each flow in isolation does not dynamo kinematically. 

Another reason to extend the studies to time-dependent flows is that 
steady flows cannot account for the irregularity of geomagnetic reversals. 
The induction equation is linear with eigenfunctions b that change with t 
only in magnitude, when 9(o") = 0, or oscillatory solutions that reverse 
polarity with fixed period 2n/S$(a). Geoma gnetic-type reversals require 
changes in the flow. ISarson and Jones! jl999) described irregular reversals 



that occurred in simulations with their 2^-dimensional model. The mecha- 
nism could be interpreted kinematically, and reversals were observed to occur 
when flu ctuations in the fl o w lead to a reduced meridional circulation. More 
recently Wicht and OlsenI (|2004h studied reversals in a fully self-consistent 



but quasi-periodic system. The reversal mechanism they proposed also ap- 
peals largely to kinematic principles and appears to reverse with approx- 
imately fixed period even when nonlinearity through the Lorentz force is 
omitted. 

The class of steady flow s explored in I— III was originally prescribed by 
Kumar an d Roberts ( 19751) and, with parameters chosen to mimic flows near 
the limit of Braginskvl (j l964). was shown to be capable of dynamo action. 



Dependence of the dynamo on a much wider range of parameter values 
was later found in I. The Kumar-Roberts flow is confined to the sphere 
of unit radius, the exterior of which is assumed to be perfectly insulating. 
Three components of the flow represent a differential rotation, a meridional 
circulation and a convective overturn, 

u = e t? + eis^ + ( e2 s| c + e3si s ). (2) 

Following the nomenclature detailed in I, the ej are constrained such that 
62 = £3 and the kinetic energy of the flow is unity. The flow is parameterised 
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by (D, M)-space, where D = D(e ), M = M(e 1 ) and \D\ + \M\ < 1. The 
parameters D and M are measures of the differential rotation and meridional 
circulation respectively. 

For a steady forcing flow, writing h(t) = e CT *B, where B is independent 
of t, Q can be expressed as the eigenvalue problem 

aB = TB. (3) 

In paper I dyn amo action was establ i shed f or approximately half the (D, M)- 
space (Fig. 12) . ISarson and Gubbinsl (^996) and III found a number of oscilla- 
tory solutions for steady flows in a region which corresponds to the dynamo 
wave solutions of the aw equations in the Braginsky limit — \D\ — » 1 in a 
manner such that 1 — \D\ 2 = c\M\ where c is a constant. The oscillatory 
region was found in I to extend only for a very narrow range in M, shown 
schematically in Fig. The majority of solutions are steady. Given the 
narrow range for M, it is apparent that only a small degree of meridional 
circulation is required to stabilise the field to steady solutions. On the other 
hand, the existence of oscillatory solutions for low M appears to be a fairly 
robust feature as the range in D for which they exist is large, and extends 
well beyond the limit of Braginsky. 

In this work, the exploration above is extended to the dynamo action of 
flows with D = D(t) and M = M(t) periodic in time, with a given period 
T. The induction equation © can be written as 

d t b = T(t)b, (4) 

with periodic forcing T(T + t) = T{t). It follows from Floquet's theorem 
(see that solutions may be written in the form b(T + t) = e aiT h(t) 
where the real part of a\ is the net growth rate over one period. Setting 
b(t) = e ffli B(t), so that B(T + t) = B(t), substitution into © defines the 
eigenvalue problem for B, 

<riB = (jr-d t )B. (5) 

The critical magnetic Reynolds number for which the field is marginally 
stable, 5R(cj) = 0, is denoted R^. 

Both the steady and non-steady eigenvalue problems © and (jHJ) per- 
mit solutions for four linearly independent spatial symmetries, axial dipole, 
axial quadrupole, equatorial dipole and equatorial quadrupole. Symmetry 
selection in the steady case was studied in II. Here only the geophysically 
interesting axial dipole symmetry will be considered. 



2 Numerical method 



Steady flows have been studie d using extensions of the method first devel- 
oped by Bullard and Gellmanl ( 1954^ . Toroidal and poloidal potentials for 
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the magnetic field are expanded in spherical harmonics, with truncation at 
degree L. A finite difference scheme is applied on N r points in the radial 
dimension leading to the discretised eigenvalue problem 

(jB = EB. (6) 

The matrix E has dimensions N r Nh x N r Nh, where after symmetry consid- 
erations the number of harmonics ~ | L 2 . As the finite difference scheme 
only connects neighbouring points, £ is block banded where each block has 
size Nh x Nh- Eigenvectors are then calculated by either by inverse itera- 
tion or by the Implicitly Restarted Arnoldi Method (IRAM) on the inverse. 
Due to the performance of both methods with respect to the distribution 
of the eigenvalues, both operate on the inverse and require the (banded) 
LU factorisation of E. Memory requirements scale like several times N 2 N ri 
depending on the stencil size of the finite difference approximation. Solu- 
tions have generally been calculated with second order differences, and L 
not much larger than twenty. The storage requirement for the large matrix 
is the limiting factor for the calculation. 

For the time-dependent eigenvalue problem © with the same spatial rep- 
resentation, applying a Fourier expansion in time introduces at least another 
factor Nt to the storage requirements. This can be minimised by permit- 
ting only sinusoidal forcings, but due to the structure of the matrix memory 
requirements are prohibitive with respect to calculation of the LU factorisa- 
tion (a few times N 2 N 2 Nt). Storage is a significant difficulty in multiplying 
by the inverse or in calculating the inverse of a suitable preconditioner for 
the time-dependent problem. 

Instead we have adopted a method that does not require storage of the 
matrix, which we call the matrix-free Krylov subspace method. It is an 
adaptation o f a method used to fi nd steady solutions of the Navier-Stokes 



"aP (Il994h Per iodicity 
( Verhulstl . ll99(ih . Writi 



equations by Ed wards et alJ (1 1994'). Periodicity of the flow is incorporated 



in the following manner ( Verhulsti . |1996|) . Writing the discrete form of (j3J) 



as 

d t b = F(t)b, (7) 

the matrix G(t) satisfying dtG{t) = F(t) G(t), with G(0) = /, is the funda- 
mental matrix of the system (J7J). Evolution of a starting solution is then 
given by 

b(t) = G(t)b(0). (8) 

For any T-periodic F(t), there exist matrices P(t) and E such that the fun- 
damental matrix can be written 

G(t) = P(t)e Et , (9) 

where E is independent of t and P(t) is T-periodic (Floquet's theorem). It 
follows immediately that the change in the solution over one period is given 
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by 

b(T) = G(T) b(0), G(T) = e tT = -4. (10) 

The stability of solutions to ® is determined by the eigenvalues A of the 
constant matrix A. If b(0) is an eigenvector of A with eigenvalue A = e CTlT , 
we find that b(T + t) = e aiT b(t) for any t. The real part of the Floquet 
exponents o\ correspond to growth rates of the solutions. Although A is 
unknown, from (jSJ) we see that the effect of multiplying by A is equivalent to 
the result of timestepping through one period. Therefore we do not have to 
calculate and store A explicitly. Note that for a steady forcing F, the period 
T can be chosen arbitrarily. 

The eigenvalue problem for A is well suited to the Arnoldi process ( ArnoldiL 



1951), which tends to pick out eigenvalues isolated in the complex plane. 
The many decaying modes have A clustered about the origin, marginal 
modes correspond to |A| close to unity. At each iteration we add to the 
Krylov-subspace given by spanjb, Ah, A k ~ l b} which we hope contains 
our solutions. In exact arithmetic the fc th Krylov subspace is equivalent to 
spanjbi, bfc} where the basis vectors b^ are obtained from the Arnoldi 
method. Numerically the latter set is better suited to span the space. The 
Arnoldi process is summarised as follows: (1) Take a suitable normalised ini- 
tial basis vector bi = lo / 1 1 1> 1 1 2 - (2) At the k th iteration evaluate (timestep) 
bfc+i = Ahk- (3) The result bfc + i is then orthogonalised against previous 
vectors in the basis by the modified Gram-Schmidt method: 

hjk = (bfe+i, bj), b k+1 := b k+1 - h jk bj] j < k. 

(4) Setting h k+ i k = ||bfc + i H2, the process continues from (2) with the next 
basis vector b k +i = bk+i/hk+i,k- Construct Q k = [bi,...,b k ] and H k = 
[hjm]j m <k- From steps (3) and (4) we expect hk+i jk — > 0. In this case, 
looking carefully at the steps above, the results of the Arnoldi process can 
be written AQ k = QkHk- Multiplying on the right by eigenvectors x of H k 
we find that they are related to those of A by b = Qjtx. For non-zero hp.+ik 
eigenvectors have residual \\Ab — Ab||2 = |/ifc+i.fc| \xk\i where x k is the last 
element of the fc-vector x. Thus, at each iteration eigenvalues A of H k are 
approximate eigenvalues of A. 

In practice the residual |/tfc+i k\\ x k\ tends to overestimate the error, and 
in our calculations the Arnoldi process is stopped when the largest eigen- 
values of H k are sufficiently converged. The number of iterations required 
is typically of order 100 or less, and so the eigenvalues of the small matrix 
/-/fc can be cheaply calculated by the QR algorithm. The memory required 
to store the basis vectors scales like NhN r k. It is possible to restart the 
Arnoldi process with a more suitable starting vector obtained from the pro- 
cess so far, but without completely restarting the new process from scratch. 
This implicit restarting allows further reduction of m emory requirements b y 
reducing the number of basis vectors at each restart ( Lehouca et al. . 1998). 
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With k small, restarting was not found to be necessary, however. The stencil 
of the finite difference scheme does not explicitly affect the storage require- 
ments. Basis vectors were therefore timestepped with a fourth order finite 
difference scheme. Timesteppi ng was performed with th e benchmarked code 
of Gibbons, Jones and Zhang ((Christensen et all 

Another advantage of the matrix- free method is that, given a timestep- 
ping code, only a few extra lines of code are required to incorporate the 
Arnoldi process, therefore leaving significantly less room for error. The 
code was verified by comparison with the matrix method used in I— III for 
the steady problem, adjustments for the periodic case in our matrix-free 
method are minimal. Table ^ a) shows the leading two eigenvalues for the 
of the steady Kumar-Roberts flow, (D,M) = (0.98354915,0.0001632689), 
at R m = 1000 calculated using the method in I— III. Table ^ b) shows the 
same eigenvalues calculated using our method. The higher order radial dif- 
ferences used in the timestepping code leads to more rapid convergence with 
N r . Table^c) shows that incorporating the Arnoldi method accelerates con- 
vergence relative to timestepping alone (and calculation of more than one 
eigenvalue is possible). The same starting vector was used for both calcula- 
tions. The period T can be chosen arbitrarily for the steady flow case, but 
if chosen too small more iterations are required and therefore more basis 
vectors must be stored. For these calculations we set T = 0.001. 

As the structure of the eigenfunctions varies with D and M, so does 
the convergence with N r and L. For most of the following analysis a radial 
resolution of N r = 50 and a spectral truncation of L = 18 was adopted. 
Checks at higher resolutions were calculated where growth rates were found 
to vary rapidly with the parameters. 

3 Results 

Periodic flows are defined by a closed orbit in (D, M)-space. We restrict 
ourselves to simple sinusoidal variations in time with a single frequency lo: 



It is the aim of this section to assess how the amplitude of variations Ad,m 
and frequency u affect the growth rates and therefore the dynamo action. 

3.1 Magnetic growth rates for time- varying flows 

Figure [21 shows growth rates for different amplitudes of variations about 
the point (Dq,Mq) = (0.25,-0.14), marked A in Fig. which lies on a 
line of minimum R m extending from the Braginsky limit point (1, 0) (see I, 
Table 5). The majority of neighbouring points have lower growth rates for 



D(t) 
M{t) 



Dq + Ad cos(wt), 
M + A M s\n(ujt). 



(11) 
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the given value of R m = 87. Figure [2] shows that the effect of oscillations of 
the flow on the growth rates is more pronounced with increased oscillation 
amplitude. 

For a steady flow, given any small real number e > there exists a finite 
time t such that transients are reduced so that |er — <j{t)\ < e, where a is the 
growth rate corresponding to the steady flow at each point on the orbit, and 
a(t) is the growth rate of an arbitrary initial field as it adjusts to the flow. 
Provided that growth rates are piecewise continuous (for example D(t) could 
be discontinuous but periodic, see Backus 1958), a frequency uj can always 
be selected low enough such that net growth rate over the cycle is close to 
the mean a of those on the orbit. The limit o\(uj) — > a as ui — > is observed 
in the numerical calculations. What is more interesting, however, is that 
with finite uj the dynamo can do much better than this mean, <J\(uj) > a, 
as seen in most figures for the growth rate. Figure |21 shows that the effect 
increases with R m and that the peak occurs at a frequency uj that increases 
in proportion to R m . 

Meridional sections of the magnetic field for this flow are plotted in Fig. 
0J The two times correspond closely to the points on the orbit which have the 
maximum (upper row) and minimum growth rates for steady flows (uj — > 0). 
The structure of these eigenfunctions is similar; regions of are generally 
well separated in space. The dissipation for these fields is larger for the lower 
panel where fields of opposite sign are squeezed towards the equator. For 
non-zero uj the location of the flux changes over the cycle, and at \ogu> = 2.6 
the field represents a smoothed version of the two eigenfunctions. Fewer 
small-scale features are present and the flow performs well as a dynamo 
(see Fig. At R m = 150 the peak frequency for o\ is logw « 2.6; 

taking T ~ ($/d) 2 as an approximate timescale for diffusion, if 5 ~ d/8 is an 
approximate length scale for the small scale features of the eigenfunctions, we 
find that the timescale for diffusion and for the peak flow oscillation coincide. 
The magnetic field is then smoothed effectively. Above this frequency the 
growth rate decreases again as the field is unable to respond to rapid changes 
in the flow. Spatial smoothing is lost and the field is close to steady — plots 
at the two times for logo; = 2.9 are almost identical. The field responds as 
though to a steady flow, retaining the stronger (smaller-scale) features from 
each eigenf unction. 

Figure El shows growth rates for an orbit about the point (0.5,-0.15), 
marked B in Fig. ^ which lies on the lower boundary of the region of suc- 
cessful steady dynamos. D remains constant and M varies to carry the flow 
outside the dynamo region. The time-dependent flow produces a positive 
effect on the growth rate. The spatial structure of the eigenvectors on this 
orbit is similar to that of the previous point considered, with well separated 
regions of positive and negative azimuthal field. For R m = 150 there is 
dynamo action only for frequencies logo; ~ 2.6, and the average growth rate 
around the orbit is negative. 



8 



We now describe a case where meridional circulation is greater than 
differential rotation, (D,M) = (-0.10,-0.45), marked C in Fig. UJ This 
point is close to where the critical magnetic Reynolds number for steady 
flows is at a global minimum. Once again time dependence of the flow 
enhances the growth rate, Fig. El A small rise can be seen in the growth rate, 
although less significant relative to the increase associated with an increase 
in the magnetic Reynolds number. Being close to the point where is 
a minimum, the fields are already relatively large scale and the smoothing 
effect of oscillations therefore has a smaller effect on the growth rate. 

Growth of the time-dependent solutions is not always found to be better 
than the mean of the growth rates on the orbit. Figure [7] shows the case 
(D,M) = (0.50,0.11), marked D in Fig. ^ where oscillations in the flow are 
initially damaging to the dynamo. If the oscillations are sufficiently rapid, 
however, the dynamo is again able to perform better. Meridional sections are 
shown in Fig. |HJ The regions of strongest flux are located very close together 
on the equator. Over the cycle radial shifts of the clover-leaf pattern of flux 
causes considerable overlap of opposite signs. This occurs mostly towards 
the outer edge of the equatorial region. As the flow oscillates, in the lower 
plot for logo; = 2.13 it can be seen that there is cancellation of flux in the 
outer region. Some flux remains at the other regions where the signs for the 
two eigenfunctions do correlate. This cancellation of fluxes over the cycle 
leads to reduced growth rates. When the oscillation is much faster, however, 
the dynamo does not have time to generate flux of opposing sign. The field 
is more steady for logo; = 2.60 and has a larger growth rate. 

The radial field in the above has little structure of interest. It is con- 
centrated mainly on the axis with opposite sign in each hemisphere (see I, 
Fig. 9a, b). The structure does not change appreciably over the cycles. In 
strength, however, it is observed to wax and wane. 



3.2 Dynamo wave solutions 

Merid ional circulation has been seen to play a key part in reversals. Sarson and Jonesl 
(1999) have studied a system in which irregular reversals are linked to a 
drop in meridional circu lation, leading to a preference for oscillatory fields. 
Wicht and Olser] (|2004l ) have recently studied a reversal mechanism that 



involves an advection of reversed flux by a large-scale SI flow. Reversals oc- 
cur quasi-periodically in their model. This behaviour may be related to the 
dynamo wave solutions obtained in III. The oscillation has the form of a dy- 
namo wave in which flux migrates along the longitudes defined by the down- 
wellings of the convective parts of the flow, which could partially explain 
the observed tendency for virtual geomagnetic poles to trac k around the 
Pacific during polarity transition ( Gubbins and SarsonL 1994). The steady 
flow model can only oscillate periodically, but we can construct a more geo- 
physically realistic reversal by choosing a time-dependent flow that traces 
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an orbit in (D, M)-space that strays into the dynamo wave region, depicted 
by the line E in Fig. ^ for a fraction / of its period. The field behaviour 
will depend on the frequency of the dynamo wave, lod and the time spent 
by the flow in the oscillatory regime. If uj 3> fuj£> the flow will only spend a 
brief time in the oscillatory regime and we expect only a minor change in the 
magnetic field. If uj <C fujjj the flow spends a long time in the oscillatory 
regime and we expect the field to oscillate several times before becoming 
steady again. The interesting case is when uj « fuio, when the field may 
only have time to oscillate for one or a half cycle, producing an excursion 
or a reversal respectively. 

We now explore reversal behaviour using periodic flows. The structure 
of the eigenf unctions for steady flows changes appreciably as M crosses zero. 
It was seen in the previous section that, where this is the case, fluctuations 
are not necessarily good for the dynamo. Instead, an orbit is chosen to 
enter the oscillatory range from the negative side. Consider the flow de- 
fined by fixed D = 0.7 and M varying sinusoidally between —0.0140 and 
—0.0020 (E in Fig. This orbit spends approximately one third of the 
time within the band of oscillatory solutions reported in III, which lies be- 
tween M = —0.0057 and —0.0010. The dynamo wave frequency for steady 
flows increases with R m and appears to saturate at about 9(<r) = 17 (see 
III, Figs 2,3); it appears to be limited by the diffusion time. Here, uj for 
the time-dependent flow must be chosen comparable with this frequency to 
give a single reversal, a value which is too low to assist the dynamo action 
significantly. 

Figure |H1 shows the the growth rate for the time dependent flows as 
a function of frequency uj. At R m = 700, 9(<r) is approximately 10 for 
steady flows in the oscillatory range. Reversing solutions may be expected 
for Tf > 7t/Sj(o"), or equivalently u < 20/3, logo; < 0.8 as the time in the 
oscillatory range is approximately one third of the cycle. Reversing solutions 
(dashed curve) are observed for uj larger than the dynamo wave frequency, 
although growth rates fall quickly when the period of the flow oscillation 
is too short to be compatible with the period of the oscillatory solution. If 
the time within the oscillatory region leads to only a half-complete reversal, 
the field exiting the region bears little resemblance to the entering field, 
which is much like the eigenfunction for these low u, and therefore leads 
to reduced growth rates. For greater uj the field does not spend sufficient 
time within the oscillatory region to reverse (solid curve) and at higher uj 
the growth rates are increased. The magnetic energy for a typical reversing 
solution is plotted in Fig. El showing a drop as the solution passes through 
the oscillatory region where the reversal occurs. It is possible to vary M so 
that / is less than a third. However, for the reversing solution in Fig. 1101 
the smooth growth rate curve, while outside the oscillatory region, indicates 
the field quickly becomes independent of the period within the oscillatory 
region, apart from in sign, due to the slow period of the flow oscillation. 
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The reversal sequence for B r at the surface is shown in Fig. 1111 Patches 
of reversed flux appear at low latitudes, strengthen and migrate polewards 
replacing th e flux at high latitudes. Th e reversal looks very similar to that 
reported by Gubbins and S arson! ( 1994h who found the pole paths during 
the reversal correlate well with the longitudes of these flux patches, located 
180° apart. If the frequency of the oscillation is too high the field is simply 
disrupted by the short period in the oscillatory region, as seen in the energy 
in Fig. 1101 This may lead to the type behaviour seen in geomagnetic excur- 
sions. Figure ^] shows that reversed patches emerge but have insufficient 
time to migrate polewards before dissipating. They still still weaken the 
dipole, however. 



4 Conclusions 

We have devised a new matrix-free Krylov subspace approach to solving the 
time-dependent stability problem that is most effective in exploring kine- 
matic dynamo action of periodic flows. It is computationally efficient, uses 
far less storage than conventional methods, and requires rather little new 
coding once time-step and eigenvalue routines are available. 

Time variation of the flow can sometimes, but not always, enhance dy- 
namo action. At low frequency the growth rate of the time dependent flow 
approaches the average growth rate for the steady flows along the orbit. At 
moderate frequency the time dependent flow can smooth out any concentra- 
tions of magnetic flux generated by the component steady flows. This can 
produce enhanced dynamo action (higher growth rate than the average) if 
the flux concentrations are isolated and of one sign. Dynamo action is possi- 
ble at certain frequencies even when the average growth rate for steady flows 
around the cycle is negative and the cycle contains mainly steady flows that 
do not generate magnetic field. The growth rate appears to be capped by 
the highest growth rate of any steady flow on the cycle. At high frequency 
the magnetic field does not have time to adjust to time changes in the flow 
and becomes almost stationary. 

Time variation does not always enhance dynamo action. When the gen- 
erated field has flux concentrations of different signs close together, the 
smoothing effect tends to destroy flux. The dynamo enh ancement for thes e 
large scale flows is not as dramatic as those reported bv iGog et al. (1999). 
Flows with similar eigenfunctions tend to have similar growth rates at the 
same R m , limiting the effect of the time variation. Where the eigenfunction 
changes dramatically dynamo action is usually impaired unless the period 
of the flow is short compared with the diffusion time. 

Solutions have been found that reverse when M is low and are associ- 
ated with the steady-flow oscillatory s olutions found in I I . This result is in 



common with the reversals studied by Sarson and Jones! (|l999l ). which oc 



11 



cur irregularly due to a drop in meridional circulation. Although flows in 
their calculations are also predominantly equatorial antisymmetric, a small 
but increased flow across the equator is observed during a reversal. It is 
unclear that this results in sufficient advection of flux to influence their re- 
versal mechanism. A large circulation exterio r to the tangent c ylind er is 
required in the reversal mechanism studied by Wicht and Olserl ((20041 ) . It 
is needed to transport reversed flux originating from plumes that protrude 
the tangent cylinder. However, it is difficult to decipher what part fluc- 
tuations would play in their model, and in particular to what degree this 
would affect the quasi-regularity of their reversals. We have shown that for 
a reversal to occur the drop in M must persist long enough for the field 
to reverse, which for this class of flows this is approximately a tenth of a 
diffusion time, or approximately the dipole decay time. A significant drop 
in magnetic energy is observed during the reversal. This arises because of 
the change in eigenfunction between the steady and oscillatory modes of the 
steady solutions. 
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M o.o 




Figure 1: Regions of dynamo solutions for steady flows and axial-dipolc 
symmetry for the flow parameterisation: D, differential rotation; M merid- 
ional circulation. The thick line running into the right hand apex contains 
oscillatory solutions related to the dynamo waves of the Braginsky limit. 
Other lines and circles describe periodic flows mentioned in the text. 
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Table 1: Comparison of computed growth rates for the K-R flow, 
(0.9834915,0.0001632689), at R m = 1000; (a) leading two eigenvalues com- 
puted using the matrix and its LU factorisation; (b) eigenvalues calcu- 
lated by the matrix- free method; (c) comparison with simple timestepping, 
T = 0.001, t k = kT and N r = 75, L = 16. 
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Figure 2: Growth rates o\ vs. u> for different amplitudes Ad = Am = A 
about the point (0.25,-0.14). 
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Figure 3: As Fig. EJfor various R m ; A D = A M = 0.1, A in Fig. □ 
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uj — > log (j = 2.6 logo; = 2.9 




Figure 4: Meridional sections, B^, = 7r/2; i? m = 150, (D,M) = 
(0.25,-0.14), Ad = Am = 0.1. Eigenfunctions were calculated for steady 
flows at the appropriate point and are plotted with independent contour val- 
ues; for uj 0, contours are plotted for the same values at different times. 
At \ogio = 2.6 the field structure looks smoothed. At logo; = 2.9 the flow 
changes too quickly for the field to respond and the field is almost steady. 
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Figure 5: Growth rates a\ vs. to for an orbit about the point (0.5, —0.15); 
A D = 0, A M = 0.1 (Fig.d B). 
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Figure 6: Rising growth rates for an orbit about (—0.1, —0.45); Ad = Am = 
0.1 (Fig. HJ C). 
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Figure 7: Growth rates for increasing uj about the point (0.5, 0.11); Ad = 0, 
A M = 0.1 (Fig.m D). 
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lo^O logo; = 2.13 logo; = 2.60 




Figure 8: Meridional sections, 4> = 0; R m = 150, (D, M) = (0.5,0.11), 
Ad = 0, Am = 0.1. The times t = ^ T and y| T correspond to maximum 
and minimum M(t) respectively. There is a radial shift of the 'clover' pat- 
tern near the equator for the eigenfunctions (u — > 0). Closely proximity of 
opposing flux leads to cancellation seen towards the outer boundary, lower 
panel with logo; = 2.13. At logo; = 2.60 the field is more steady. 
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Figure 9: »(er), R m = 700 at (D, M) = (0.7, -0.008), A D = 0, A M = 0.006 
(Fig. [IJ E). For the dashed curve 9(cri) = u/2. 
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Figure 10: Magnetic energy (minus net growth) for parameters as Fig. 03 
Reversing solution logcj = 1.30 (solid); failed reversal logo; = 2.00 (dashed). 
Vertical bars represent the period in the oscillatory regime. 
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Figure 11: B r at the surface during the reversal in Fig. I1UI Patches of 
reversed flux near the equator migrate polewards, replacing the flux at high 
latitudes with reversed field. 




Figure 12: B r at the surface during the failed reversal in Fig. 1101 Reversed 
flux patches have insufficient time to migrate polewards before dissipating. 
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